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Abstract 

Based on the color dipole QCD approach we discuss the multiple scattering 
mechanism (MSM) of suppression of back-to-back azimuthal correlations of 
high px forward dijets in DIS off nuclei. We quantify the effect in terms of 
acoplanarity/decorrclation momentum of jets. For hard dijets the decorrelation 
momentum is found to be of the order of the nuclear saturation momentum 
Qa- Minijets with the transverse momentum below the saturation scale are 
proved to be completely decorrelated. The analysis of DIS indicates that in 
heavy ion collisions the contribution of MSM to the reactions dynamics can be 
of the order of magnitude of the decorrelation effect observed at RHIC. 
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1 Introduction 



Recent observation of gradual disappearance of azimuthal back-to-back corre- 
lations of high pt particles with centrality of Au Au collisions at RHIC 1) is 
presently viewed as a consequence of the jet energy loss in hot quark-gluon 
plasma produced in central collisions. Another possible explanation is that 
initial and final state interactions associated with multiple parton scatterings 
on hard stage of the process produce uncorrelated monojets. 

To get an idea of the size of the multiple scatterings (MS) effect at RHIC 
we start with the QCD description of the breakup of photons into forward dijets 
in small- a: deep inelastic scattering (DIS) off nuclei in the saturation regime 

m. In ® we reported a derivation of the general formula for the two-body 
transverse momentum distribution. Our formalism, based on the technique 
El El EJ)^ mc i uc i e s consistently the diffractive attenuation of color dipoles and 
effects of transitions between different color states of the gg-pair propagating 
through the nucleus. Opacity of nuclei brings in a new scale Qa which separates 
the regimes of opaque nuclei and weak attenuation El El El EJ p or hard 
dijets the decorrelation momentum is found to be of the order of the nuclear 
saturation momentum Qa ^* ■ For parton momenta below the saturation scale 
Qa the evolution of sea from gluons was shown to be dominated by the anti- 
collinear, anti-DGLAP splitting EJ. As a result, minijets with the transverse 
momentum below the saturation scale are proved to be completely decorrelated 

El) 

Turning back to the RHIC observation we argue that MS mechanism may 
contribute substantially to the effect of vanishing of back-to-back azimuthal 
correlations of high px hadrons in central nucleus-nucleus collisions. 



2 Breakup of photons into hard dijets on nuclear targets 

We consider DIS at x < xa = I/Rah^n 1 which is dominated by interactions 
of qq Fock states of the photon and make use of the conventional approximation 
of two t-channel gluons in DIS off free nucleons EH> , The two-gluon exchange 
approximation amounts to neglecting diffractive DIS off free nucleons which 
is justified by a small fraction of diffractive DIS, rjn < 1 HJ. We sum 
unitarity cuts of the forward Compton scattering amplitude which describe the 
transition from the color-neutral qq dipole to the color-octet qq pair. It should 
be emphasized that we are interested in the unitarity cuts which correspond to 
the genuine inelastic DIS with color excitation of the target nucleus. 

For x < xa the propagation of the qq pair inside nucleus can be treated 
in the straight-path approximation. Let b + and b be the impact parameters 
of the quark and antiquark, respectively, and SU(b-|_,b_) be the S-matrix for 



interaction of the qq pair with the nucleus. Regarding the color states |c) = 
|1), |8) of the qq pair, we sum over all octet and singlet states. Then the 2-jet 
inclusive spectrum is calculated in terms of the 2-body density matrix as 

, n dam n = -i-- / d 2 b'd 2 W rf 2 b+d 2 b_ 
dzd 2 p+d 2 p_ (2tt) 4 J + - + 

x exp[-zp + (b+ - b+) - ip (b - b'.)]^** 

A* c 

-(1;A\s a \a-i)(1;A\s a \A;1)}. (i) 

In $ = fy(Q 2 , z,h + — b_) stands for the wave function of the qq Fock 
state of the photon with virtuality Q 2 and the photon light-cone momentum 
fraction z carried by the quark. Notice, that the calculation of the 2-body 
density matrix enter four straight-path trajectories S^(h±) , <5( 2 )(b' ± ) and Sa 
and S A describe the propagation of two quark-antiquarks pairs, qq and q'q' , 
inside a nucleus. In the integrand of Q we subtracted the coherent diffractive 
component of the final state. 

Upon the application of closure to sum over nuclear final states A* the 
integrand of can be considered as an intranuclear evolution operator for the 

2-body density matrix (for the related discussion see refs. EJ^EH) 
^^< J 4|{<l|51|c)}|A*)<A*|{<c|^|l)}|A) = 

A* c 

= {A\l^2(l\S* A \c}{c\S A \l)^\A). (2) 

Let the eikonal for the quark-nucleon and antiquark-nucleon QCD gluon ex- 
change interaction be T£xO) and T°x(b), where T£ and T± are the SU(N C ) 
generators for the quark and antiquarks states, respectively. The vertex V a for 
excitation of the nucleon g a N — > N* into color octet state is so normalized that 
after application of closure the vertex g a g b NN is S a t. Then, to the two-gluon 
exchange approximation, the S'-matrix of the (gg)-nucleon interaction equals 

S N (b + , b_) = 1 + z[r;x(b + ) + Tl X (b-)]V a - \[Tlx(h+) + ^ X (b_)] 2 . (3) 

The profile function for interaction of the qq dipole with a nucleon is T(b + , b_ ) = 
1 — 5Ar(b + , b_). For a color-singlet dipole (T? +T") 2 = and the dipole cross 
section for interaction of the color-singlet qq dipole with the nucleon equals 

a(b+-b_) = 2 / d 2 b + (7V|r(b+,b_)|iv) 



JV?-1 



|d 2 b + Wb + )-x(b_)] 2 . (4) 



2iV c 

The nuclear S-matrix of the straight-path approximation is 

A 

SU(b+,b_) = l[S N (b + b,,b -bj), 

3 = 1 

where the ordering along the longitudinal path is understood. We evaluate the 
nuclear expectation value in (J2J in the standard dilute gas approximation. To 
the two-gluon exchange approximation, per each and every nucleon Nj only 
the terms quadratic in x(bj) must be kept in the evaluation of the single- 
nucleon matrix element (Nj\S^SN\Nj) which enters the calculation of S a Sa- 

Following the technique developed in El 0* we can reduce the calculation of 
the evolution operator for the 2-body density matrix J5J) to the evaluation of 
the S'-matrix S4A(b+, b_, b' + , b'_) for the scattering of a fictitious 4-parton 
state composed of the two quark-antiquark pairs in the overall color-singlet 
state. Because (T?)* = — T", within the two-gluon exchange approximation the 
quarks entering the complex-conjugate S A in can be viewed as antiquarks, 
so that 

J2(MS*a\c)(c\S a \1) = J2^c(c'c\S 4A \U) , (5) 

c c'c 

where S^a (b' + , b'_ , b + , b_) is an S- matrix for the propagation of the two quark- 
antiquark pairs in the overall singlet state. While the first qq pair is formed by 
the initial quark q and antiquark q at impact parameters b + and b , respec- 
tively, in the second q'q' pair the quark q' propagates at an impact parameter 
b'_ and the antiquark q' at an impact parameter b^_ . 

If ai(b' + , b'_, b + , b_) is the color-dipole cross section operator for the 4- 
body state, then the evaluation of the nuclear expectation value for a dilute gas 
nucleus in the standard approximation of neglecting the size of color dipoles 

compared to a radius of heavy nucleus gives 

S 4A (b' + , b'_ , b+, b_) = cx P { - ^<74(b' + , bl, b+ , b_)T(b)} , (6) 

where T(b) = J db z nA(b z , b) is the optical thickness of a nucleus at an impact 
parameter b = j(b + +b' + + b_ +b'_) , and n A (b z ,b) is nuclear matter density 
with the normalization J <i 2 bT(b) = A . The single-nucleon S'-matrix © con- 
tains transitions from the color-singlet to the both color-singlet and color-octet 
qq pairs. However, only color-singlet operators contribute to (Nj\SjfSN\Nj), 
and hence the matrix 04 only includes transitions between the |11) and 
color-singlet 4-parton states, the 1 18} states are not allowed. 



After some color algebra, we find matrix elements an = { 1 1 1 cr^ 1 1 1 ) , dig = 
(T8i = (11 |cT4 188) and ergs = (88|cr 4 |88) (details of the derivation are presented 

in EJ). For forward hard jets with the momenta 3> Q\, which are produced 
from interactions with the target nucleus of small color dipoles in the incident 
photon the two eigenvalues of the operator 04 are 



and 



£2 ~ cm 



For small color-singlet dipoles £2 ~ cm ~ and the nuclear distortion factor 
takes on a simple form 

^^(^|(l|^|c)(c|^|l)|A)-(l;^|^|A;l)(l;A|^|^;l) 

A' c 

= ((11| + VA^(88|) exp |-ia 4 T(b)| |11) - exp j-i«7 U r(b) 

» y/Nf^T— (l-exp 



(7) 



The introduction of this nuclear distortion factor into Q gives the hard dijet 
inclusive cross section: 



dai. 



d 2 bdzd 2 p+d 2 A 



nb) g« (b ,,-)/ d v/o> (i - K)5i ^ 



(8) 



where the acoplanarity/decorrelation momentum A is 

A = p+ + p , 

the probability of finding j spatially overlapping nucleus in a Lorentz-contracted 
nucleus is 

^ (b;j) _l7(, + l,2Wb)) 



J- 



2\ c v A {h) 



1 



v A {h) = -a s a T(b), 
\ c = N 2 J{N 2 C - 1) 

and x) = L dyyi^ 1 exp(— y) is an incomplete gamma- function. The func- 
tion 

/ (i) («)= /nd a « i /(«i)*(«-E K< )' / (o) («)=%) w 



in eq. JEJ is a collective gluon field of j overlapping nucleons introduced in LL2» . 

The probabilistic form of a convolution of the differential cross section on 
a free nucleon target with the manifestly positive defined distribution /w (k) in 
|(HJ| can be understood as follows. Hard jets originate from small color dipoles. 
It is the quantum mechanical interference that suppresses interaction with soft 
gluons of the small-sized color-singlet qq state. However, the first inelastic 
interaction inside a nucleus converts the qq pair into the color-octet state, in 
which color charges of the quark and antiquark do not neutralize each other, 
rescatterings of the quark and antiquark in the collective color field of intranu- 
clear nucleons become uncorrelated, and the broadening of the momentum 
distribution with nuclear thickness follows a probabilistic picture. 

3 Decorrelation of dijets in DIS off nuclei: numerical estimates 

To quantify the azimuthal decorrelation of two forward jets we find it convenient 
to introduce the mean transverse acoplanarity momentum squared (Aj.(b)), 
where At is transverse to an axis of the jet with higher momentum. It is 
assumed that jets are hard, |p+| ^> Qa- The convolution property of the hard 
dijet cross section (JSJ suggests 



broadening term The sign sa in l|l(J|) reflects the kinematical limitations on 
p and k in the practical evaluation of the acoplanarity distribution. In a 
typical final state it is the harder jet with larger transverse momentum which 
defines the jet axis and the acoplanarity momentum A will be defined in terms 
of components of the momentum of softer jet with respect to that axis, for 
instance, see . For the sake of definiteness, we present numerical estimates 
for the Gedanken experiment in which we classify the event as a dijet if the 
quark and antiquark are produced in different hemispheres, i.e., if the azimuthal 
angle (j> between two jets exceeds ir/2, the quark jet has fixed |p + | and the 
antiquark jet has higher transverse momentum |p + | < |p_| < 10|p+| (in the 
discussion of the experimental data one often refers to the higher momentum 
jet as the trigger jet and the softer jet as the away jet ^'). The free-nucleon 
quantity (A t )jv can be estimated starting with the small- A expansion for 
excitation of hard, p+ ^> e 2 = z(l — z)Q 2 , light flavor dijets from transverse 
photons 




daw 



-e 2 f a em a s (p 2 + ) [z 2 + (1 - zf] 

7T 

I dG(x,A 2 ) A 2 



dzd 2 p + d 2 A 



x 



A 4 ' dlogA 2 (e 2 + p+)(e 2 + p+ + A 2 ) ' 



(11) 
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Figure 1: The mean acoplanarity momentum squared (A r ) n for DIS on a free 
nucleon target with production of trigger jets with the transverse momentum 
higher than p+ for several values of Q 2 . The numerical results are for x = 0.01 
and the input unintegrated gluon SF of the proton is taken from ref. [3. 



where dG{x, A 2 )/31og A 2 = T{x, A 2 ) is the unintegrated gluon structure 
function (SF) of the free nucleon HU> , Then a quick estimate 



(A 2 T ) N 



2 H*,P+) 

P+ 2G(x,p%)> 



(12) 



correctly describes the numerical results shown in fig. 1. As far as the dijets 
are hard, p 2 ^ > z(l — z)Q 2 ~ \Q 2 , the acoplanarity momentum distribution 
would not depend on Q 2 , which holds still better if one considers o~t + o~l- 
This point is illustrated in fig. 1, where we show (A 2 ^) at at z = 1/2 for several 
values of Q 2 . Because of this weak dependence on Q 2 here-below we make no 
distinction between DIS and real photoproduction, Q 2 = 0. In the practical 
evaluations of the nuclear contribution (^(b))^ one can use the eq.© which 
gives the result 



(K 2 T (b)) A « -A c Qi(b) 



log- 



OTcQA(b) 



(13) 



where 



12 /kn _ 



Q\{h) » ^-a s (Q 2 )G(z,Q 2 )T(b) . 



(14) 
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Figure 2: The atomic mass number dependence of nuclear broadening contri- 
bution. (K^(h)) a, to the mean acoplanarity momentum squared for real pho- 
toproduction off nuclei at x = 0.01. The input unintegrated gluon SF of the 
proton is taken from ref. 



For average DIS on heavy nuclei the reference value is (Q\ u (b)) = 0.9 

(GeV/c) 2 , see EJ. The atomic mass number dependence of nuclear broadening 
(k^)a for jets with p + = 4 GeV/c in average DIS off nucleus is shown in fig. 2. 
The principal reason why (k^)a is numerically small compared to (Q\ u {b)} 
is that even for such a heavy nucleus as 197 Au the no-broadening probability 
in average DIS is large, (wAu(b, 0)) ~ 0.5. A comparison of the free nucleon 
broadening (A^)jv from fig. 1 with the nuclear contribution (/^(b))^ from 
fig. 2 shows that the nuclear mass number dependence of azimuthal decorrela- 
tion of dijets in average DIS off nuclei will be relatively weak. 

However, nuclear broadening will be substantially stronger for a subsam- 
ple of central DIS events at b ~ 0. In fig. 3 we show for the gold, 197 Au, target 
a dependence of the averaged nuclear broadening (nj,(b))A on the impact pa- 
rameter at several values of p+ . 

One can enhance Q\ and nuclear contribution (/^(b))^ still further se- 
lecting DIS events when the photon breaks up into the qq pair on the front face 
of a nucleus. Experimentally, precisely such events are isolated by selecting 
very large multiplicity or very high transverse energy of produced secondary 
particles ( and references therein) . 

Now we comment on the recent finding by the STAR collaboration of 
a disappearance of back-to-back high pr hadron correlation when going from 

peripheral to central gold-gold collisions at RHIC . The application of 
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Figure 3: The impact parameter dependence of the nuclear broadening contribu- 
tion, (Kj,(b))^. to the mean acoplanarity momentum squared from peripheral 
DIS at large impact parameter to central DIS at b = for several values of the 
away jet momentum p + . The numerical results are for x = 0.01 and the input 
unintegrated gluon SF of the proton is taken from ref. 



color dipole formalism to hard hadron-nucleus interactions » suggests that 
our analysis of acoplanarity of forward hard jets can be readily generalized to 
mid-rapidity jets. One only has to choose an appropriate system of dipolcs, 
for instance, the open heavy flavor production can be treated in terms of the 
intranuclear propagation of the gluon-quark-antiquark system in the overall 
color-singlet state. At RHIC energies jets with moderately large pr are for the 
most part due to gluon-gluon collisions. In our language that can be treated 
as a breakup of gluons into dijets and azimuthal decorrelation of hard jets 
must be discussed in terms of intranuclear propagation of color-octet gluon- 
gluon dipoles. For such gluon-gluon dipoles the relevant saturation scale Q\ A 

is larger 1^ than that for the quark-antiquark dipoles by the factor 2A C = 
Ca/Cf = 9/4. Arguably, in central nucleus-nucleus collisions distortions in the 
target and projectile nuclei add up and the effective thickness of nuclear matter 
is about twice of that in DIS. Then, the results shown in fig. 2 suggest that for 
central gold-gold collisions the nuclear broadening of gluon-gluon dijets could 
be quite substantial, (/cf.(b = 0)}auAu ~(3-4) (GeV/c) 2 for average central 
AuAu collisions and even twice larger if collisions take place at front surface of 
colliding nuclei. 

The principal effect of nuclear broadening is a reduction of the probability 



of observing the back-to-back jets 



(A 2 



PQ,) « - "~ T ' N „ (15) 
V ' {n T {h)) A + (A 2 T ) N 1 ' 

and one needs to compare (A t )n to («y(b))yi . Our eq. (|12|l for the free nucleon 
case holds as well for the gluon-gluon collisions. Then the results shown in fig. 1 
entail (A t )n ~ (k t (0)}auAu ~(3-4) (GeV/c) 2 at the jet momentum p + = 
pj = (6-8) GeV/c and our nuclear broadening would become substantial for all 
jets with p + below the decorrelation threshold momentum pj. In practice, the 
STAR collaboration studied the azimuthal correlation of two high-p^ hadrons 
and for the quantitative correspondence between the STAR observable and 
azimuthal decorrelation in the parent dijet one needs to model fragmentation 

of jets into hadrons (for the modern fragmentation schemes see ^»), here 
we notice that the cutoff p + in our Gedanken experiment is related to the 
momentum cutoff pT,m%n of associated tracks from the away jet, whereas our 
jet of momentum p_ can be regarded as a counterpart of the trigger jet of 
STAR. The STAR cutoff px = 2GeV/c corresponds to the parents jets with the 
transverse momentum p + ~ (2 — 3)px =(4-6) GeV which is comparable to, or 
even smaller than, the decorrelation threshold momentum pj — (6-8) GeV/c. 
Then eq. I|15fl suggests that in the kinematics of STAR the probability to 
observe the back-to-back away and trigger jets decreases approximately twofold 
from peripheral to central AuAu collisions, 

P(0) « 0.5, 

and perhaps even stronger, so that our azimuthal decorrelation may contribute 
substantially to the STAR effect. 



4 Remark on breakup of photons into semihard dijets 

In our previous analysis of the single particle spectrum it has been discov- 
ered that the transverse momentum distribution of sea quarks is dominated 
by the anticollinear, anti-DGLAP splitting of gluons into sea, when the trans- 
verse momentum of the parent gluons is larger than the momentum of the sea 
quarks. That suggests strongly a complete azimuthal decorrelation of forward 
minijets with the transverse momenta below the saturation scale, p± < Qa- 

In GJ> this limiting case has been considered in detail. The principal point is 
that the minijet-minijet inclusive cross section depends on neither the mini- 
jet nor decorrelation momentum. This observation proves a disappearance of 
the azimuthal correlation of minijets with the transverse momentum below the 
saturation scale. 
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